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1. Introduction 



Phase transitions in various spin models have served as wonderful examples for the study of 
quark-hadron transition in quantum chromodynamics (QCD) and the related SU (N c ) gauge theo- 
ries, where N c is the number of colours. While average magnetization serves as the order parameter 
in the former case, the Polyakov loop, L, defined as the product of the timelike gauge links at a given 
site, is the order parameter for the deconfinement transition [jl]|. On an Euclidean Ng x N T lattice 
L(x) is defined at a site x as 

L(x) = ^TrU^ =l U 4 (x,x Q ), (1.1) 

where U^(x) are the gauge variables associated with the directed links in the jitfh direction, pt = 
1,4. As in the spin models again, it is convenient to define its average over the spatial volume, 
L = Y*xL{x) /N a . (1^1) was use d to establish a second order deconfinement transition in numerical 
simulations of the SU (2) pure gauge theory. Since then it has been used for similar studies of 
the deconfinement phase transitions for a variety of N c [Q], for establishing the universality [^] 
of the continuum limit, as well as for theories with dynamical quarks [ffl|. Further, the predicted 
universality [||] of critical indices has also been numerically verified [Bj. Indeed, one hopes to be 
able to construct effective actions |0] for L in a Wilsonian RG approach. These will be similar to 
the spin models in the same universality class but with possibly additional interaction terms. A 
large number of models of quark-hadron transitions use the Polyakov loop as the order parameter 
for the deconfinement transition as well. 

An order parameter should be physical, i.e., independent of the the lattice size. This is indeed 
so for spin models for sufficiently large lattices. For SU (N c ) gauge theories, this requirement 
means in addition independence from the lattice spacing a in the continuum limit. Furthermore, it 
must be so in both the phases it seeks to distinguish. As is the case for any bare Wilson loop, the 
Polyakov loop, needs to be renormalized for this to be true. Since the bare Polyakov loop is further 
known to decrease progressively with N t , suggesting it to be zero in the continuum limit in the high 
temperature phase, renormalized L is even more desirable to have. 



2. Results 



The physical interpretation of the order parameter as a measure of the free energy of a single 
quark, (L(T)) = exp(— Fq{T)/T) provides a straightforward clue for renormalization. Since many 
years various attempts to remove the divergent contribution in the single quark free in the contin- 
uum limit have been made. These include computations employing lattice perturbation theory [g|, 
use of the heavy quark-antiquark free energy [[J, fits to (L) on A^ T -grids JTo| ] and an iterative direct 
renormalization procedure [ 11 ] for (L) among others. 

Here I advocate Jl2| ] another, perhaps better, method to define renormalized (L). Let me 
elaborate why this maybe so. The definition of Ref. ^ needs heavy quark potential at short 
distances. Lattice artifacts are at their worse when one is at such short distances, with maximal 
violation of the rotational invariance. Finite volume of the lattice also enters in defining the large 
distance between the heavy quarks, or Polyakov loops. Similarly the iterative procedure used in 
Ref. [11] to obtain the renormalization constants needs large lattices in both spatial and temporal 
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Figure 1: The average Polyakov loop as a function ofT/T c for four different scales. The lattice sizes are as 
indicated in the key. 

directions. Physically perhaps an undesirable aspect of the definition of Ref. [|ll]] is that it works 
only on the plasma side, i.e., for T > T c , where T c is the position of the peak in the Polyakov 
loop susceptibility. The definition has so far been employed only in the T >T C for pure gauge 
theories for which L is an order parameter. It would clearly be nice if the renormalization procedure 
is applicable to the usually employed (\L\), which is used as an order parameter on finite volumes. 

I obtain a renormalized Polyakov loop which is valid for both the phases below and above T c 
Jl2|]. It can be defined in any spatial volume, and it becomes the true order parameter in the infinite 
spatial volume limit. Of course, it is also physical, i.e., /^-independent on finite volumes as well. 
Indeed, it seems to work rather well for a range of temporal lattice sizes, including N t > 4. I use the 



fixed scale approach [13] to do so. It was introduced to minimize the computational costs for the 



zero temperature simulations needed to subtract the vacuum contribution in thermodynamic quanti- 



ties such as the pressure and to isolate pure thermal effects in computation of T c []14p. Furthermore, 
its advantage is that all the simulations stay on the line of constant physics in a straightforward 
way. What I argue is that it is indeed this advantage which also permits an easy renormalization 
of the Polyakov loop. Although these considerations are general, and apply to any SU (N) gauge 
theory as well as any quark representation, I shall consider below the simplest case of the SU (2) 
lattice gauge theory to illustrate how and why it works. 

Recall that the temperature T is varied in this approach by varying N t , holding the lattice 
spacing a, or equivalently the gauge coupling j8 = 2N c /g fixed. The single quark free energy 
F[,(N T ,a) is then obtained from the L by the canonical relation, 

ln(\L\) = -aN z F b {N r ,a) . (2.1) 

The subscript b reminds us that one obtains the bare free energy this way. If the chosen coupling is 
j8 c , corresponding to the position of the peak of the |L|-susceptibility in the usual fixed N T approach, 
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Figure 2: The heavy quark free energy F a function ofT/T c for two different scales. The lattice sizes are as 
indicated in the key. 

and it lies in the scaling region, then the physical deconfinement temperature T c = 1 /N x ^ c a c , and 
T /T c = N ZjC /N T in the fixed scale approach, with the free energy given by Fb(T /T c ,a c ). Writing it as 
a sum of a divergent and a regular contribution, one has a c Fb(T /T c ,a c ) = a c F(T /T c ,a c ) —a c A(a c ), 
where A is the divergent free energy in physical units. Clearly, the divergent contribution will be 
same at all temperatures in the fixed scale approach since it depends only on a c . 

Since j8 c , or a c , is known precisely for the Wilson action of the SU (2) theory for many different 
N T , I chose four different scales labelled T c $, T c ^, T c %, T c \2 corresponding to the known transition 
couplings on N x = 4, 6 JB|] and 8, 12 [||] respectively: j3 cl = 2.2991, fi c2 = 2.4265, j3 c3 = 2.5104, 



and j8 c 4 = 2.6355. Note that T/T c is given simply by n/N T with n = 4, 6, 8 and 12 respectively. 
Employing then Af T = 3 to 12, 1 varied the temperature in the range 2>T/T C > 0.6. Note that fixed 
scale a c leads to a constant spatial volume in physical units in each case in contrast to the usual 
fixed N t approach where the spatial volume varies with T. I used a variety of spatial lattice sizes. 

Figure [j] shows the results for the thermal expectation value of L as a function of the temper- 
ature in the units of T c . In most cases, I used both a random and an ordered start. The errors are 
corrected for autocorrelations. The agreement in the data for the two starts suggest the statistics of 
200K iterations to be sufficient. As expected, the four different scales, T C 4 , T c ^ , T c $ , and T c \2 lead 
to four different curves for the order parameter. One also sees the known feature of (L) — > as 
a c — > even in the deconfined region. Figure displays the behaviour of the bare free energy for 



just two scales, obtained by using the eq.([2.1|). The stars are for the scale T c g while the squares are 
for r c 6 corresponding to the higher lattice spacing of the two. The figure reinforces the expectation 
of the effect of the divergent free energy, since the free energy increases with the decrease in the 
lattice cut-off a c . 

Any two different scales, a c \ and a C 2 have their respective divergent contributions, a c \A(a c \) 



and a C 2A{a C 2). Multiplying eq.(2. 1) by Afy, for j = \ and 2 corresponding to the critical N x for the 
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Figure 3: Heavy quark free energy F as a function ofT/T c with a constant shift, as explained in the text. 
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Figure 4: Renormalized Polyakov loop versus T /T c using the shifted free energy of the upper figure, as 
explained in the text. 

scale choices above, i.e, 6 and 8, one obtains 

T Hm = _w^i t (22) 

where Fb(T /T c ,a c j) /T c = F{T /T c ,a c j)/T c —A(a c j)/T c . Thus the free energies at the same temper- 
atures but two different scales are related by a mere constant, [A(a c i) — A(a c z)]/T c . For the four 
scales considered here, this implies 3 such constants. Figure |3] shows the results for the free energy 
with three constant shifts in the free energy determined by demanding coincidence at the highest 
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Figure 5: Renormalized Polyakov loop versus T /T c after subtracting the would-be divergent contribution. 

T = 2T C . A universal curve for the free energy seems to result as a result for a wide range of 
T > T c . The results for the low temperature phase are seen to be volume dependent, as expected. In 
the infinite volume limit, the free energy should increase to infinity in the confined phase whereas 
it should essentially remain constant in the deconfined phase. Such an expectation is indeed borne 
out by the results in the Figure [3| For the same physical volume, the free energy appears to be 
a-independent in the T <T C phase as well, as seen by comparing the crosses and the stars. 

Finally, it should now be clear how one can obtain a universal curve for the order parameter 
from the universal free energy curve. The (\L\) corresponding to scale j8 C 2 should simply be multi- 
plied by the factor exp(./V T [A(a c i)/r c — A{a c 2)/T c ]) and then the date will lie on a universal curve. 
This is exhibited in Figure ^ for all the four scales. It is worth noting that the same universal order 
parameter results in both below and above T c by fixing only three constants for the four scales ex- 
hibited. The entire low and high temperature region of the order parameter is uniquely fixed, and 
appears to be universal. 

From the Figure Q it appears as though the approach of (\L\) to unity is slow and from below. 
It is, however, known since long [ |I7| ] that perturbation theory predicts L — > 1 from above at very 
large T : L = 1 +C^g 3 + 0(g 4 ), where ct,(N c ) > is a constant. The solution to this apparent 
paradox can be traced to the usual fact that a renormalized quantity depends on the scale chosen 
to define the scheme for renormalization. In my case, the inclusion of a constant free energy 
A(a c )/T c for the chosen scale a c defines the choice. The details of the shape of the physical order 
parameter are therefore scale-dependent in the plasma phase but it is universal none the less once 
a choice is made. Moreover, any further change of scale leads to a computable change in the 
shape. Indeed, in order to mimic the perturbative renormalization scheme I estimated the point- 
divergent contribution. At the highest temperature 2T C , I fitted the results for the four scales to 
— ln(|Ly|) = F(2T C )/2T C +B -N X j/2. Having thus determined the coefficient B of the would-be 
divergent contribution at the scale T c n, I eliminated the B-dependent contribution at that scale. The 
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renormalized L at the other three scales were related to it by the same shifts as before. Figure || 
displays the resultant (|L|). It crosses unity at about 1.5 T c . Since all the multiplying factors tend 
to unity at large T, the approach of this L to unity is from above at large T. Note also that large 
spatial volumes, aspect ratio of ~ 10, are needed for this L ~ in the low T phase. 

3. Summary 

In conclusion, I showed that the fixed scale approach leads to a natural definition of a physical, 
^-independent, order parameter which is defined in both the confined and the deconfined phases. 
The definition itself does not depend on any lattice artifacts or the lattice size in the deconfined 
phase , and works very well for even coarse lattices (a < l/4T c ). Moreover, it displays the ex- 
pected behaviour in the confined phase as the physical volume is increased, suggesting that the so 
determined physical free energy of a single quark in the confined phase, F, goes to infinity in the 
infinite volume limit. Eliminating the point-divergent contribution leads to a high temperature be- 
haviour consistent with perturbation theory. It is straightforward to generalize this idea to SU (N c ) 
gauge theories and QCD as well as to sources in higher representations. 
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